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ABSTRACT 


*An  Z  -invertible  nontinite  totally  positive  matrix  A  is  shown  to  have 


one  and  only  one 


^ain 


diagonal’’.  This  means  that  exactly  one  diagonal  of  A 


has  the  property  that  all  finite  sections  of  A  principal  with  respect  to 


this  diagonal  are  invertible  and  their  inverses  converge  boundedly  and 
-1 


entrywise  to  A  .  This  is  shown  to  imply  restrictions  on  the  possible  shapes 
of  such  a  matrix.  In  the  proof,  such  a  matrix  is  also  shown  to  have  an 
invertible  LUU  factorization.  In  addition,  decay  of  the  entries  of  such  a 
matrix  away  from  the  main  diagonal  is  demonstrated.  It  is  also  shown  that  a 
bounded  sign-regular  matrix  carrying  some  bounded  sequence  to  a  uniformly 
alternating  sequence  must  have  all  its  columns  in  Cq  . 
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SIGNIFICANCE  AND  EXPLANATION 

Spline  approximation  is  often  most  effective  when  the  breakpoint  (knot) 
sequence  can  be  chosen  suitably  non-uniform.  At  the  same  time,  the  standard 
approximation  schemes  (such  as  least  squares  approximation,  or  interpolation 
at  suitable  interpolation  points  by  splines)  are  so  far  only  known  to  be 
usable  and  bounded  as  long  as  the  breakpoint  sequence  is  almost  uniform.  The 
problem  of  showing  existence  and  uniqueness  of  bounded  spline  approximants  to 
bounded  data  boils  down  to  showing  invertibility  of  a  certain  infinite  matrix 
A.  The  distinguished  feature  of  this  matrix  is  its  total  positivity,  i.e., 
all  minors  of  A  are  nonnegative.  In  this  paper  we  show  that  an  arbitrary 
invertible  totally  positive  matrix  has  a  main  diagonal  and  one  only.  By  this 
we  mean  that  there  is  exactly  one  band  or  diagonal  with  the  property  that 
finite  sections  of  the  matrix  whose  main  diagonal  is  a  piece  of  that  band  are 
invertible  and  their  inverses  converge  to  the  inverse  of  the  matrix  as  the 
section  becomes  large.  This  says  that  it  is  possible  to  approximate  the 
solution  of  the  infinite  system  by  the  solution  of  appropriately  chosen  finite 
sub-systems. 

These  results  illustrate  the  power  of  total  positivity,  since  such 
results  fail  to  hold  for  arbitrary  infinite  matrices.  It  is  hoped  that  the 
better  understanding  of  total  positivity  of  nonfinite  matrices  obtained  in  the 
present  report  and  its  predecessors  will  aid  in  the  resolution  of  certain 
outstanding  problems  concerning  the  stability  of  spline  approximation  schemes. 


The  responsibility  for  the  wording  and  views  expressed  in  this  descriptive 
summary  lies  with  MRC,  and  not  with  the  authors  of  this  re^xjrt. 


STRUCTURE  OF  INVERTIBLE  (BI) INFINITE  TOTALLY  POSITIVE  MATRICES 
C.  de  Boor*,  Rong-qing  Jia**  and  A.  Pinkus*** 

Introduction.  An  Investigation  into  the  structure  of  invertible  (bi)infinite  totally 
positive  (or,  tp)  matrices  was  begun  in  (B^l  and  [CDMS^],  motivated  by  certain  questions 
[B^l  [Ml  in  spline  approximation.  First  results  concerned  banded  matrices  [B2]  [CDMS2I  and 
strictly  banded  block  Toeplitz  matrices  (COHS^ ]  .  Somewhat  surprisingly,  such  restriction 
to  banded  matrices  was  found  to  be  unnecessary  in  [BFP],  where  it  was  proved  that  a  tp  (or, 
more  generally,  a  sign  regular)  matrix  A  mapping  *,„(J)  to  t„(I)  is  onto  if  and  only 
if  it  has  a  uniformly  alternating  sequence  in  its  range.  This  latter  condition  means  that, 
for  some  bounded  x  ,  sup^[Ax) (i )(ax) (i+l )  <  0 

A  similar  result  had  been  proved  in  (B^l  for  band  matrices,  following  up  a  conjecture 
in  (M).  The  argument  in  (B^l  provided  much  additional  information  about  the  inverse  of  a 
banded  tp  matrix.  In  the  present  paper,  we  obtain  the  same  information  for  an  arbitrary 
4^-invertible  tp  matrix. 

I  ^  J 

We  prove  (in  Section  1 )  that  an  f^-invertible  tp  matrix  A  6  R  has  a  "main 
diagonal",  by  this  we  mean  that  all  sufficiently  large  square  submatrices  of  A  taken 
from  consecutive  rows  and  columns  and  principal  with  respect  to  a  fixed  diagonal  are 
invertible  and  their  inverses  are  bounded  uniformly.  This  implies  that  A"’  is  the 

pointwise  limit  of  these  inverses.  We  even  show  that  such  convergence  is  monotone.  As  a 

J  —I  T  K  K^K 

consequence,  u'^A”  or  its  negative  is  again  tp,  with  b  6  R  the  diagonal  matrix 

tor  which  oNk)  =  (-1)*^  ,  all  k  e  K  . 

In  addition,  we  show  (in  Section  2)  the  uniqueness  of  such  a  main  diagonal.  For  the 
proof,  we  establish  a  result  of  independent  interest,  namely  the  existence  and  uniqueness 
of  an  invertible  LbU  factorization  tor  A  .  By  this  we  mean  the  possibility  of  writing  A 
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as  the  product  LUU  ot  bounded  and  boundedly  invertible  tp  matrices,  with  L  and  its 
inverse  unit  lower  triangular,  D  and  its  inverse  diagonal,  and  U  and  its  inverse  unit 
upper  triangular. 

We  also  follow  up  some  consequences  of  the  existence  of  a  main  diagonal  for  the 
matrix  A  :  (i)  A  can  have  only  certain  shapes.  Explicitly,  we  must  have  J  «  I  +  r  for 
some  r  .  We  settle  (in  Section  3)  the  related  question  of  possible  choices  for  I  and 
J  in  case  we  only  know  that  A  is  tp  and  onto,  (ii)  The  main  diagonal  must  be  bounded 
away  from  zero.  This  leads  to  the  observation  (made  in  Section  4)  that  the  entries  of  A 
must  decay  at  least  linearly  away  from  the  main  diagonal.  We  leave  open  the  intriguing 
question  whether,  in  fact,  this  decay  is  exponential,  as  was  shown  for  arbitrary  matrices 
with  a  banded  Inverse  in  (D]. 

In  Che  last  section,  we  settle  a  conjecture  from  [BFPl  whose  consideration  gave  much 

•  IM 

impetus  to  the  present  paper.  We  prove  that  a  sign  regular  f^-invertible  matrix  A  6  R 
must  already  carry  Cq(J)  to  Cq(1)  ,  i.e.,  all  its  columns  must  vanish  at  Infinity. 

We  use  the  following  notation  and  conventions. 

We  use  lower  case  letters  to  denote  elements  of  R^  ,i.e.,  real  valued  functions  (or, 
sequences)  on  some  integer  interval  (or,  more  generally,  some  integer  set)  I  ,  with 
x(i)  the  ich  entry,  or  value  at  1  ,  of  the  sequence  x  .  If  K  is  a  subset  of  I  , 
then  x,^  denotes  the  restriction  of  x  to  K  .  The  i-th  unit  sequence  is  denoted  by  e^ 
and  defined  by 

e^(j)  6^^  ,  all  i,j. 

Further,  tor  x  S  R^  , 

|xl(i)  :»  |x(i)|  ,  all  i  . 

We  use  Che  customary  abbreviation  4.(1)  for  the  normed  linear  space  of  all  sequences  x 
on  1  wi  th 

Ixl  :»  sup  |x(i)| 

i 

finite,  and  use  c^d)  for  its  closed  linear  subspace  ot  all  bounded  sequences  x 
vanishing  at  infinity,  i.e.,  for  which  {i  «  I  :  |x(i))  >  a}  is  finite  for  all  positive 


a  .  v»e  use  t^{I)  for  the  norned  linear  space  of  all  absolutely  summable  sequences  x 
on  I  ,  i .  e . ,  wi  th 

Ixl^  :»  S  lx(i)| 

1 

finite,  and  we  make  use  of  the  fact  that  *„<I)  is  the  dual  of  i^(I)  ,  which,  in  turn, 

is  Che  dual  of  Cq(I)  . 

Let  also  J  be  an  integer  interval  (or,  more  generally,  a  set  of  integers).  We  are 
Ix  J 

interested  in  A  6  R  which  carries  i^(J)  to  4^(1)  ,  and  we  Identify  such  a  matrix 

A  with  the  linear  map  from  4„(J)  to  4^(1)  induced  by  it.  Such  A  is  necessarily 
bounded  and  we  denote  by  lAl  its  map  norm.  If  A  is  one-one  and  onto,  it  is  boundedly 
invertible.  We  will  then  say  that  A  is  1^-invertible,  for  short.  As  with  sequences, 

|A|{i,J)  :=  lA(i,j)|  ,  all  i,j. 

We  write  A(*,j)  for  the  sequence  which  makes  up  the  j-th  column  of  A  ,  and  A(i,  •)  for 
the  i-th  row.  We  call  the  sequence  AC'.'-frl  the  r-th  band  or  diagonal  of  A  . 

If  KXL  ^  ixj  ,  then  we  denote  by  Aj^  the  restriction  of  A  to  KxL  .  At  times, 
it  IS  also  convenient  to  use  the  alternate  notation 


and  even 


A(K]  :=  A^  . 

On  replacing  square  brackets  by  round  brackets,  we  get  to  the  determinant: 

A^  ...  J  :*  detA^  ...  j  . 

Thus  A{i,i)  and  A(i)  both  denote  the  i-th  diagonal  entry  of  A  ,  -  the  latter  involving 
an  abuse  of  notation  since  A({i})  would  have  been  correct. 

We  recall  that  A  is  sr^  (:»  sign  regular  of  order  k)  if,  tor  some 
€^,  ...,  W  {-I,l}  and  tor  all  KxL  ^  I*J  with  |k|  =  |l|  «  k  , 

where  |k|  denotes  the  cardinality  of  the  index  set  K  .  If  all  1  ,  then  we  use  tpj^ 

(:=  totally  positive  of  order  k)  instead  of  st)^  .  It  A  is  stj^  tor  all  k  ,  then 

A  IS  simply  called  sr  (tp).  We  intend  to  make  use  of  Sylvester's  determinant  identity 
(see,  e.q.,  [Kj  p.3|)  and  Madamard’s  inequality  (see,  e.g.,  (K;  P.U81). 
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1-  Existence  of  a  o^in  diagonal*  v^e  establish  that  any  ^^-invertible  natrix  A 
must  nave  a  main  «liaijonai*  Hy  tnis  we  mean  that  A  is  tne  stable  limit  ot  its  tinite 
sections  wnich  are  principal  witn  resj>ect  to  a  particular  diagonal,  it  is  this  diagon<ai 
that  we  then  call  the  "main  diagonal”.  We  use  here  the  definite  article  in  anticipation  of 
the  uniqueness  established  in  the  next  section, 

ixj 

^eorea  1*  If  the  matrix  A  6  R  is  tp  and  l_»invertibie,  tnen,  tor  some  r  ,  and 

for  all  finite  intervals  M  ^  1  ,  A|^  is  invertible  and  (A^j  converges , 

monotonely  in  each  entry,  to  A"'  a^  M  — 1  , 

Proof*  since  A  is  onto,  there  exists  a  bounded  sequence  x  such  that 

<Ax)<i)  =  all  1  . 

With  this#  the  argument  for  Theorem  1  of  (BFP]  and,  especially,  the  Claim  proved  there, 
provide  us,  for  each  finite  interval  L  c  i  ,  with  an  index  set  K  J  so  that 
*A,  fc-  '  *  ^  •  (In  the  Claim,  choose  S  =*  A  ,  hence  S*  ,  choose  I  =  h  , 

choose  u(i)  =  (-1)^.  ail  i,  and  choose  for  s  an  extremal  for 

F^:  Za^A(i,*)  t— -e  Ea^v(i)  with  v  =  u  .) 

„  -1  L  Jxi 

Extend  each  fo  a  matrix  C  fe'  R  by  taking  its  value  to  be  zero  oft 

k*L  .  Since,  tor  each  1  ,  IC^(‘,i)l  <  Ixl  ,  we  are  permitted  to  choose  a  sequence  ot 

integer  intervals  L  approaching  I  for  which,  for  each  i  I  ,  converges 

*  i 

weak  (i.e.,  pointwise  on  )  to  some  y  €  ,  Then,  tor  each  m  «  I  ,  the 

sequence  A(n,*)  oeinq  in  &^(J)  ,  we  have 

li  L  fe‘(m),  meb'l 

l,Ay  Jim)  =  Urn  Z  A(m,s)C  (s,i)  =  lim  (  \  =  e  (m) 

L*1  L*I  I  0  ,  f4lI 

✓ 

snowing  that  Ay^  =  e^  ,  therefore  =  a“\*,i)  ,  all  i  .  We  conclude 

converges  pointwise  on  and  so,  in  particular  pointwise,  to  a”^(*,i)  . 

The  remainder  of  the  proof  relies  on  the  following  known  fact. 
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Lena  1.  if  a  v  n 


IS  tp  dnd  xnvertibie,  then,  tor  any  interval  M  <=  {l 


n)  , 


so  is  C  :=  B[M]  =  Bji  ,  and 

0  <  (-)'^*^c’’(i,j)  <  (i,j)  .all  i,j  e  M  . 


The  nonsinqularity  of  C  is  a  consequence  of  Hadamard's  inequality.  The  rest  of  the 
lemma  folloirfs  by  repeated  application  of  the  special  case  M  =  {l,...,n-l}  .  For  this 

case,  use  the  checkerboard  nature  of  the  inverse  of  a  tp  matrix  and  the  rank-one  correction 
formula  for  an  inverse  (called  the  Woodbury-Snerman-Morrison  formula  by  Numerical  Analysts) 
which  in  this  case  states  that 

C"’(i,j)  =  B"'(i,j)  -  b"’ (i,n)B"’ (n, j )/B~^ (n)  ,  i, j=1 , . . . ,n-1 . 

Of  course,  this  identity  also  follows  from  Cramer's  rule  and  Sylvester's  determinant 
identity  using  pivot  block. 

As  a  corollary,  we  find  that,  tor  an  invertible  tp  B  6  r"*"  ,  1/B(i)  <  B~^(i)  ,  or, 

1  <  B(i)B  '(i)  ,  all  i  ,  hence 

<  B(i)  ,  all  i  . 

This  shows  that  the  (main)  diagonal  entries  of  an  invertible  tp  matrix  are  bounded  away 
from  zero. 

Returning  now  to  the  invertible  matrices  Aj^  obtained  earlier,  we  conclude  that 
their  (main)  diagonal  entries  must  all  be  bounded  below  by  V*xl  .  Assuming  for  simplicity 
of  notation  that  0  «  I  ,  it  follows  that,  tor  every  L  containing  0  ,  the  index  for 

which  A(0,j^)  lies  on  the  (main)  diagonal  of  Aj^  must  come  from  the  set 
{j  S  J  :  A(o,j)  *  1/lxl)  ,  and  this  set  is  finite  since  A(0,  •)  ®  .Vie  are  therefore 
permitted  to  choose  a  further  subsequence  of  intervals  L  increasing  to  I  (and  all 
containing  0  )  tor  which  ]^  =  r  for  some  fixed  integer  r  .  Since  j  and  L  are 
intervals  and  is  square,  this  implies  that  L+r  ^  J  for  all  such  L  . 

Since  ,  or,  more  precisely,  its  extension  to  all  of  Jxi  converges 

-1  L 

pointwise  to  A  yet  C  (],*)  =0  tor  all  j  jJ  K  ,  it  follows  that  every  j  «  J  is 
eventually  in  every  K  .  Thus,  for  every  interval  M  ^  I  ,  there  exists  L  so  that  M+r  £ 

K  .  By  Lemma  1,  ^^.M+r  invertible  and  lA^  M+r"'*  '  *\  k"' '  ^  *’'"  '  ''bother,  by  that 
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lemma,  M  M  implies  that 

We  conclude  that  ,  or,  more  precisely,  its  extension  to  all  of  J*I  ,  converges, 

iJXI 

monotonely  in  each  entry,  to  a  bounded  matrix  C  S  K  which,  by  a  repeat  of  the  earlier 
argument,  is  necessarily  A~^  .  | | | 

Corollary  1 .  The  main  diagonal  of  A  is  bounded  below  by  lA  ^ I  ^  .  More  precisely, 
A(i,i+r )A”' (i+r ,i)  *  1  for  all  i  . 

Corollary  2.  The  matrix  (  (-)^'*'^A  '(i+r,j)  )  is  tp. 

The  second  corollary  improves  on  Proposition  1  of  (BFP)  in  which  the  same  conclusion 
rs  only  reached  under  the  additional  assumption  that  the  columns  of  A  are  all  in  Cq(I)  . 
(See  also  Section  S. ) 

Theorem  1  implies  that  ^  t_-invertible  tp  matrix  A  has  a  main  diagonal  in  the  sense 
introduced  in  tB^l  :  For  some  r  and  for  all  sufficiently  large  intervals  M  ,  is 

Invertible  and  lim  sup„.-lA„  •  <  •  .  This  means  that  one  can  approximate  the 

1  M I M+r 

solution  X  of  the  nonfinite  linear  system  Ax  =  b  by  the  solution  x  :=«  b,^  of 

the  finite  system  »  b^  ,  at  least  when  b  and  x  are  in  Cq  .  We  show  in 

Section  b  that  x  W  Cq  if  and  only  if  b  6  Cq  . 

Finally,  Theorem  1  implies  that  I  +•  r  »  J  for  some  r  .  Hence  it  is  no  restriction 

to  assume,  after  a  shift,  that  I  =  J  .  If  I  has  a  first  or  a  last  element,  this  then 

forces  band  U  to  be  the  main  diagonal. 
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2.  Existence  of  stable  LU  factorization  and  uniqueness  of  main  diagonal.  The  LDU 

tactorization  ot  a  matrix  A  of  order  n  is  a  basic  tool  in  the  solution  of  linear 
algebraic  systems.  The  matrix  A  is  written,  if  possible,  in  the  form  LDU  ,  with  L 
lower  triangular,  D  diagonal,  and  U  upper  triangular.  Both  L  and  U  are  required  to 
be  unit  triangular,  i.e.,  to  have  all  their  diagonal  entries  equal  to  1  ,  and  this  insures 
uniqueness  of  the  factorization  (see  the  proof  ot  Lemma  2  below).  Such  a  factorization  can 
be  obtained  in  case  all  upper  left  principal  minors  of  A  are  nonzero.  By  Hadamard's 
inequality,  this  condition  is  satisfied  for  an  invertible  tp  A  .  Further,  if  LDU  =  A  , 
then  (Du](i,*)  is  the  particular  linear  combination  ot  the  first  i  rows  of  A  which 
vanishes  at  1,2,...,i-1  and  in  which  A(i,*)  has  coefficient  1  .  Therefore 
(Du](i,*)  =  a(i,...,i-1,  \  ]/A( 1 , . . . ,1-1 ) 

since  the  right  hand  side  has  the  same  properties.  Also, 

A(1,...,i)  =  L(1,...,i)D(1,...,i)U(1,...,i)  =  D(1)...D(i)  , 

hence 

D(i)  =  A(1, ...,1)^(1,. ..,1-1)  . 

Therefore,  finally, 

LC,j)  =  A(l,...,j-1,  )/A(1,...,j)  . 

We  have  recalled  these  well  known  facts  to  facilitate  discussion  of  an  LDU 

factorization  of  a  non finite  matrix*  We  intend  to  obtain  such  a  factorization  by  a  limiting 

process  and  need  therefore  to  consider  the  behavior  of  the  LDU  factorization  of  A  as  we 

add  rows  and  columns  to  A  « 
ixj 

Let  now  A  €  R  .We  say  that  the  matrices  L»  D,  and  u  provide  a  (normalized) 

LOU  factorization  for  A  provided  A  =  UDU  *  L  is  unit  lower  triangular  with  band  0  as 
its  rightmost  (nontrivial)  band,  U  has  band  0  as  its  only  nontrivial  band,  and  U  Is 
unit  upper  triangular*  We  call  such  a  factorization  (boundedly)  invertible  it  each  factor 
IS  bounded  and  boundedly  invertible,  with  L**'  and  U~^  being  again  triangular. 

We  are  not  certain  that  having  a  "good"  diagonal  (i.e*,  a  diagonal  all  ot  whose 
principal  minors  are  nonzero)  is  sutficient  for  having  an  lUU  factorization;  it  may  not 
ever*  oe  sufticient  to  have  a  main  diagonal*  hut  it  A  is  also  tp,  then  having  a  main 


aiagonal  is  sutticient  for  having  a  bounuedly  invertible  LUU  factorization,  ihis  insures 
tnat  the  main  diaqonai  is  unique,  as  we  will  see. 


I«BBa  2.  _If_  L,  U,  U  and  ,  1/^  and  are  both  invertible  normalized  LOCI 

factorizations  for  A  (as  maps  from  to  ),  then  L  =  ,  D  =  ,  U  =  . 

Proof.  Let  band  r  be  the  leftmost  band  of  U  ,  and  band  r^  the  leftmost  band  of 
,  and  assume  without  loss  that  r  <  r,  .  Then,  by  assumption, 

L,~’ld  =  D,U,U"’ 

and  the  left  hand  side  is  lower  triangular  with  band  0  as  its  rightmost  band  and  equal  to 
the  interesting  diagonal  of  D  ,  while  the  right  hand  side  is  upper  triangular  with  band 
r^-r  as  its  leftmost  band,  and  equal  to  the  interesting  diagonal  of  .  Since  r,  >  r  , 
it  follows  that  r,  =  r  ,  hence  D  =  D,  ,  and  further  L^"’l  =  1  =  U^U"’  .|{| 

Theorem  2.  A  ^  t_-invertible  and  tp,  hence  has  a  main  diaqonai,  say  band  r  , 
then  A  has  an  invertible  normalized  LUO  factorization  L,  D,  U  ,  with  band  r  the 
leftmost  band  of  U  . 

Proof.  Assume  first  that  r  •  0  .  further,  assume  without  loss  that  A^  := 

Al-n,...,nl  is  a  typical  principal  section  tor  that  main  diagonal  of  A  .  Let  L^, 

Ujj  be  the  LUU  factorization  tor  A^  .  Then 

U^(i)  »  A(-n, . . . ,1 )/A( -n, . . . ,i-1 )  «  1/A(-n, . . . ,i J (1 ) 

and ,  by  Lemma  1  and  Theorem  1 , 

A(-n,...,i)"’(i)  <  '  A'’(i) 

and  A[-n, . . .  ,i  I~\i )  increases  as  n  grows.  We  conclude  that 

lu  "’l  <  lA  "’l 

n  n 

and  decreases  nonotonely  to  some  invertible  diagonal  matrix  D  for  which 

i)(i)  >  lA  '  I  '  ,  all  l  .  Therefore,  U  is  boundedly  invertible.  Niext  consider  n 

saw  earlier  in  this  section  that 


(d  U  ](i,d)  =  A(-n . i-t,  '■  ]/A(-n,...,-a-l)  . 

V  n  n'  j 


Use  Sylvester's  determinant  identity  on  A[-n-l # . < • , i-l #  ^  J  with  pivot  block 
A [ -n, . • . , i-1 ]  to  conclude  that 


a( -n-1 , . . . , 1-1 ,  ^  )a( -n, . , , ,i-l ) 


A(-n'l 


i— 1)A^— n#***/!— 1#  } 


a('"  ;•••<> 

^ -n-1 , .  • ,  1-1  ^-n  3  * 


which,  by  tp,  implies  that 

We  conclude  that  decreases  monotoneiy  to  some  matrix  V  ,  necessarily  tp  and  upper 

triangular,  theretore  converges  to  the  tp  unit  upper  triangular  matrix  U  :=  D”'v  • 

Further,  <  A^  ,  hence 

n  n  n 

■  u  I  <  lA  IlD  "’l  <  lA  Ha  "’i  ,  thus  lul  <  IaHa"'|  . 

n  n  n  n  n 

The  corresponding  argument  shows  that  converges  to  some  tp  unit  lower  triangular 

matrix  L  with  ILI  <  HaHa  '•  .  Finally,  A  ,  hence  also 

'  n  n  n  n  ' 

IL  'l  ,  lu  'l  <  lA  Ha  'i  <  IaHa”' •  .  Indeed,  |A-~' I  is  again  tp  and  ,  because 

n  n  nn  ''n'^'' 

of  the  checlierboard  nature  of  inverses  of  tp  matrices,  give 

the  UUL.  factorization  tor  l^n”^  I  '  ‘•®*<  factorization  obtained  by  starting  Gauss 
elimination  at  the  lower  right  corner  rather  than  the  upper  left,  hence  satisfy  the  same 
bounds.  For,  if  B*(i,j)  :=  b(-i,-j)  ,  ail  i,j,  i.e.,  the  order  of  both  the  rows  and  the 
columns  is  reversed,  then  tp  is  preserved  as  well  as  the  max-rowsum  norm  while,  with  H,  F, 
V  an  LUU  tactorization  tor  B  ,  the  matrices  M*  ,  K*  ,  and  V*  provide  a  UDL 
factorization  tor  b*  .  This  shows  that  also  L  and  U  (and,  of  course  [1  )  have  band  0  as 
main  diagonal,  and  are,  in  particular,  invertible. 

It  now  r  /  0  ,  then  we  obtain  such  an  LDU  factorization  (with  each  tactor  havinu  nana 
0  as  main  diagonal)  tor  ttie  matrix  AL"*^  ,  with  E  the  shift,  (Ea)(i)  ;  =  a(i+l)  ,  all  i  . 
but  then  A  =  UJ(UE*^)  does  it  .  ||| 


Corollary.  It  also  f>.iii.<  s  of  a  is  a  main  diagonal  fur  a  ,  f.en  s  =  r  . 


-g- 


Proof.  By  Theorem  2 


A  has  an  invertible  normalized  LDU  factorization  with  band  r 


as  the  leftmost  nontrivial  band  for  the  upper  triangular  factor,  as  well  as  one  with  band  s 
as  the  leftmost  nontrivial  band  of  the  upper  triangular  factor.  But,  by  bemma  2,  these  two 
factors  must  be  the  same.  | | | 

3.  The  shape  of  a  nonfinlte  tp  matrix.  The  fact  that  any  invertible  tp  matrix 

ixj 

A  U  R  has  a  main  diagonal  implies  that  I  -f  r  »  J  for  some  r  .  In  particular,  it 

rules  out  certain  combinations  of  I  and  J  which  might,  offhand,  be  possible.  For 
example,  the  choice  I  ^  -M  is  not  possible  in  conjunction  with  the  choice  J  =  N  .  If  we 
merely  know  that  A  is  onto  or  one-one,  certain  other  combinations,  though  certainly  not 
all,  are  possible. 

In  this  context,  the  only  properties  of  the  integer  intervals  I  and  J  which  matter 
are  whether  or  not  they  have  a  first  and/or  last  element.  In  effect,  there  are  just  four 
choices  for  I  and  J  : 

finite,  M,  -M,  2  . 

If  J  is  finite  then  A  cannot  be  onto  unless  |j|  >  |l|  .  The  same  conclusion  holds  in 
spirit  when  J  is  nonfinite. 

IX  J 

Theorem  3.  Let  A  «  R  be  tp  and  carry  t^(J)  onto  t„(l>  •  Let  I  and  J  each 
be  one  of  M,  -M,  or  2  .  Then  I  c  J  . 

Proof.  There  is  nothing  to  prove  unless  J  =  N  or  -R  .  Assume  without  loss  of 
generality  that  J  s  M  .  Assume  by  way  of  contradiction  that  I  >  -M  or  S  .  A  being  onto, 
we  may  pick  a  bounded  x  for  which  Ax  »  ((-1)*^)  .  Pick  any  i  e  I  .  Then,  since 
A(i,*)  6  t^(J)  ,  we  nay  pick  n  so  that 

j  >  n  implies  A(i,j)  <  Ixl”'  . 

Since  I  a  -M  or  Z  ,  the  interval  L  :=  li-n,il  is  in  I  .  By  (BPPl  (as  used  at  the 
beginning  of  the  proof  for  Theorem  1  above),  we  may  pick  a  subset  K  of  J  for  which 
lA  'is  Ixl  .  but  then,  with  K  =:  {k  <  ...  <  k  )  ,  we  must  have 


-10- 


A(i,fc  )  >  Ixl  ^  yet  k  >  n  , 
n  n 

and  this  is  a  contradiction  to  our  choice  of  n  .  |  |  | 

4a  Uecay  a%iay  froa  the  ealn  diagonal.  The  fact  that  an  l^-lnvertlble  tp  matrix  A 
has  a  main  diagonal  insures,  by  Corollary  1  of  Theorem  1,  that  it  has  a  diagonal  which  is 
bounded  away  from  zero.  Ibis  implies  at  least  linear  decay  away  from  that  diagonal.  More 
precisely,  we  have  the  following 

jxi 

Theorem  4.  Let  A  6  R  be  tp^  and  carry  1^  to  itself.  If  d  :=  inf  A(i)  >  o  , 
then  sup^  A(i,i+s)  =  0(1/|sj)  a£  js)  — *  “  . 


Proof.  Let  s  >  0  .  Since  A  is  tp2  ,we  have 

A(i,i+lt)A(i+k,i+s)  >  A(i,i4s)A(i+k,i-»k)  ,  k=1,...,s 


or 


s 

I  A(i,i+k)  max 
k»1  k«1,,.. 


min  Ad-tk) 
k=1 , . . . ,s 

This  implies  that  lAl^  ..lAl/d  >  |s| 

=  lAI^  ^  while  sup^  £)^|A(i,k)|  =  lAI 


A(i+k,i+s) 

>  S  A(l,l+S)  . 

supj^  A(i,i+s)  ,  using  the  fact  that  max^  j|A(i,})| 
.  An  analogous  argument  proves  the  inequality  for 


negative  s  .  {  |  | 


ixj 

Corollary.  _I£  A  6  R  is  tp  and  f_-invertible,  then,  for  some  r  , 
A(i,i+s)  <  lAl^  ^itCA)  /|s-r|  ,  all  i,  s  . 

Proof.  By  Corollary  1  to  Theorem  1,  we  may  choose  d  =  Ia"^ l”'  ,  hence 

lAl/d  =  IAIIa"’|  =:  K(A)  .  Ill 


In  case  a"'  is  also  banded,  we  know  from  [D)  that  the  entries  of  A  decay 
exixjnentially  away  from  the  main  diagonal  of  A  •  txplicitly  Ib^)  , 


-n- 


A( 


MA)-!  \  1  s-r  1  /m 


*  ^  -  1  ^ 


in  case  A~^  is  m-banded.  In  such  a  case,  A~^  also  decays  exponentially,  in  a  trivial 
way.  It  would  be  nice  to  know  whetner,  in  tact,  all  t^-invertible  tp  matrices  decay 
exponentially  away  from  their  main  diagonal. 


5.  Invertible  sr  aatrices.  It  is  not  possible  to  extend  the  above  results  concerning 

main  diagonals  to  invertible  sr  matrices.  This  is  illustrated  by  the  particular  matrix  B 

ZXZ 

obtained  from  the  identity  matrix  in  R  by  reversing  the  order  of  the  rows.  Every  band 
has  singular  principal  minors  of  arbitrarily  large  orders.  Of  course,  one  could  weaken  the 
notion  of  main  diagonal  to  demand  only  that  there  be  some  sequence  of  principal  submatrices 
converging  stably  to  A  .  It  is  possible  to  show  that  an  I^-invertible  sr  matrix  has  such 
main  diagonals.  For  the  particular  matrix  B  ,  though,  every  band  becomes  main  in  this 
weaker  sense,  i.e,.  the  distinction  becomes  empty.  Ibis  is,  of  course,  not  surprising 
since,  for  B  ,  there  is  really  nothing  to  distinguish  one  band  from  another. 

Some  of  the  work  reported  on  in  the  preceding  sections  was  initially  based  on  the 
realization  that  the  following  conjecture  is  true. 

Conjecture  ((BFP)).  The  columns  of  a  tp  I_-invertible  matrix  are  already  in  Cq  . 

Indeed,  this  conjecture  is  now  a  simple  consequence  of  Theorem  4.  It  is  pointed  out  in 
(BFpJ  that  the  discussion  there  of  the  inverse  of  a  sr  matrix  could  be  considerably 
shortened  if  this  conjecture  were  to  hold  for  sr  matrices.  We  now  prove  the  following 
generalization  of  the  conjecture. 

Theoreu  5>  Let  s  c_  R*^  be  a  normed  sequence  space  whicn  contains  the  unit  vector  e^ 
tor  all  j  S  J  and  contains'  |x|  it  it  contains  x  ,  and  let  A:S  — ♦  be  bounded 

and  sr.  If  Ax  uniformly  alternates  tor  some  x  6  S  ,  then  ACq  ^  Cq  ,  i.e.,  the  columns 


Proof.  We  have  to  prove  that 


K  :=  {i  ti  I  :  lA(i,k)  I  >  a} 

is  finite  for  all  k  «  J  and  all  oi  >  0  •  By  considering  -A  and/or  inverting  the  order  of 
the  rows  of  A  if  necessary,  we  can  ensure  that  A  is  tP2.  Then,  tor  j  <  k  ,  i '  <  i  ,  we 
have 

A(i,j)A(i',k)  <  A(i',3)A(i,k) 


or,  tor  1'  6  K  , 


Therefore , 


A(i,j)  < 


A(i,k) 


A(i  ' ,  j  ) 


E  |x(])|  supA(i,3)  <  £  A(i-,j)|x(3)|  <  «A(x|l 

jkk  i>i'  j<k 


or 


with 


E  |x(j)|IAC,]).„,l  <  - 

j<k  ^  ' 


(K)  :a  (inf  K,  sup  K)  . 

By  inverting  the  order  of  rows  and  of  columns  of  A  (which  preserves  tp2)»  we  obtain  the 
same  inequality  for  the  sum  over  j  >  k  ,  hence  altogether 

E  >AC,j)|^jl|x(3)|  <  -  . 

Since  3 ^  uniformly  alternates,  this  implies  that  we  can  find  some  tini te 

index  set  L  so  that  already  maps  to  a  uniformly  alternating  sequence,  and 

this,  by  a  standard  sr  result  (see,  e*g.,  (K;  Chap. 5,  ^1)  or  the  beginning  of  the  proof  of 
Theorem  1  in  (BFPJ),  implies  that  l(Kj|  <  |l|  ,  thus  |k|  <  “  .  |)| 


This  theorem  makes  it  possiDle  to  give  a  shorter  proof  of 

Ttleorea  2  of  (BFP).  A  W  is  sr  and  f_-invertihle ,  then  is  also  sr. 

by  mimicking  the  proof  of  Proposition  1  of  (UFPl,  since  the  assumptions  allow  me 
conclusion  that  A  maps  Cy(J)  to  Cy(I)  .  Hence  Corollary  2  to  Theorem  l  t't  IlKPl  shows 


with 


that  A“'  is  the  pointwlse  limit  of  certain  inatrices 
converges  to  I  .  Explicitly,  E^  equals  Aj^  on 
K  some  index  set  depending  on  L  .  This  implies  that 
limit,  a”'  • 


as  the  index 

interval 

h 

KXL 

and  vanishes 

Off 

KKL 

,  with 

eI- 

is  sr,  hence 

80  is 

its 

pointwise 

I 


I 

I 

I 
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